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Distributed Generation (DG) units are also called as Decentrahzed 



to maximize the 
selecting types, 
penetration level 



Generation and Wbedded Generation. The objective is 
DG penetration level, minimization of loss by optimally 
locations and sizes of utility owned DG units. The DG 
could be limited by harmonic distortion because of the lonlinear current 
injected by inverter-based DG units and also protection coordination 
constraints because of the variation in fault current caused by synchronous- 
based DG units. Hence the objective is to maximize DG penetration level 
from both types of DG units, taking into account power balance constraints, 
total harmonic distortion limits, and protection coordination constraints. The 
Social Learning Particle Swarm Optimization (SLPSO) algorithm is used to 
maximize the overall DG penetration level and the propose^ system is tested 
in the IEEE-30 bus system. 
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I. INTRODUCTION 

DG generally refers to small-scale (typically 1 kW - 50 
MW) electric power generators that produce electricity at a 
site close to customers or that are tied to an electric 
distribution system [1]. Optimal DG placement (ODGP) 
can improve network performance in terms of voltage 
profile, reduce flows and system losses and also improve 
power quahty and rehabihty of supply. The DG units are 
of two types namely utility owned DG units and customer 
owned DG units [2]. In the case of a utiHty owned DG 
installation, the utiHty has to optimally plan the location 
and size of the DG units in order to improve network 
benefits and reliability [3]. In the case of customer owned 
DG installation, the utility planner should conduct 
feasibility and assessment study to evaluate any tedhnicaf 
issues resulting from the new installation of customer 
owned DG installation [4]. 

A new optimization approach that employs an Artificial 
Bee Colony (ABC) algorithm to determine the optimal DG 
units size, power factor and location in order to minimize 
the total system real power loss [5]. The closed form 
equations are used for determining allowable penetration 
levels of DG resources [6]. The effect of DG on 
coordination will depends on the size, type, and placement 
of DG. The paper [7] explores that the effect of high DG 
penetration on protective device coordination. The 
advantages of installing DG in distribution systems are 
improving reliability, mitigating voltage sags, unloading 
sub transmission and transmission system [8]. 

In [9] the harmonic filter planning is presented as a 
combinational optimization problem, which is solved using 
Genetic Algorithms. Here the current injection method is 
used to determine the harmonic flow in a radial 
distribution feeder. In ref [10] a fuzzy-based approach for 
optimal placement and sizing of fixed capacity J^aill^^n: )^ t 



II. Problem Formulation 

The objective of this problem is to maximize the DG 
penetration level, by selecting the optimal size of inverter 
based DG units, using Particle Swarm Optimization (PSO) 
algorithm. The constraints of the problem include bus 
voltage limits, individual harmonic limits, total harmonic 
hmits and power balance constraints. 

A. Objective function: 

The main objective is to maximize the DG penetration 
level with respect to total system capacity. The objective 
function F can be defined as follows: 



Maximize F (%) = 



Minimize Pl = TJi^^Loss^^ 



100 



(1) 



B. Equality Constraints: 

The power balance constraints at fundamental frequency 
for each bus i can be given as follows : 
p . _|_ pifii'V _ p 
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C. Inequality Constraints: 

Total Harmonic Distortion Limits 
The total voltage harmonic distorti[)n at each bus i is as 
follows: 

ijt%)-^A^^y^-^0 < THD^^ (5) 



.CO 



(4) 



radial distribution networks in the presence of voltage and 
current harmonics. 

A new method, based on interactive Multi Objective 
Nonlinear Programming (MONLP), using Genetic 
Algorithms (GAs) to study passive filter planning [11]. 
When the system condition changes, the passive filte r s 
may be most likely to cause other power quality problem 
or have less impact on mitigation of harmonics. An 
analytical approach [12] is used to minimize the system 
losses by installing various types of DG units and finding 
optimal size and power factor of four types of distributed 
generation (DG) units. 

PSO is being used to evolve not only the network 
weights, but also the network structure. The method is so 
simple, efficient and this approach is effective for any 
network architecture [13]. The particle swarm optimizer 
has been found to be robust and fast in solving nonlinear, 
non-differentiable, multi-modal problems. 



where THD™is the maximum permissible total voltage 
harmonic distortion which is usually set at 5%. 



D. Protection Coordination: 
An inverse time over current 



characteristics is 



considered for the OCR relays as given in the following: 



m- 



(6) 



Where, 



A & B are constants that depend on the type of OCR 
characteristics. 

^^-^i is the time dial setting of OCR. 

^sc.ij is the short circuit fault current flowing through 
the relay i because of fault at location j. 

'pj is the pickup current of the relay i above which it 
starts operating. 
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III. Social Learning Particle Swarm 
Optimization 

In 1921, some British birds were first seen to open milk 
bottles in the small town of Swaythling. In the following 
25 years, such observations had been continually reported 
from numerous other sites spreading all over the Great 
Britain and even some other areas in the European 
continent. This is the first evidence of social learning, 
where the birds are believed to learn to open milk bottles 
by observations and interactions with other birds, instead 
of learning by themselves. 

During the past decades, various mechanisms have been 
proposed and discussed in social learning theory, e.g., 
stimulus enhancement and local enhancement, 
observational conditioning, contagion and 'sociaf 
facilitation. Among these mechanisms, the most interesting 
social learning mechanism is imitation, which is 
considered to be distinctive from other social learning 
mechanisms, because imitation, which operates across a 
whole community, could lead to population- level 
similarities of behavior such as culture or tradition. Such 
population-level similarities may imply convergence of a 
dynamic system, thus providing its essential applicability 
in an evolutionary algorithm [14]. 

Inspired by social learning mechanism, an imitator will 
learn the behaviors of different demonstrators in the 
following manner: 



otherwise 



(7) 



where :{t) is the j-th dimension of particle i's behavior 

vector in generation t, with t E [1,2,3,4 m} and 

; E {1, 2, 3 , 4 . . . . . . . n}; AXij it -h 1} is Sien^efiaviOT ) 

correction. Taking into account the fact that in a society, 
the motivation to learnfrom better individuals may vary 
from individual to individual (typically, better individuals 
are less willing to learn from others), we define a learning 
probability PLi for each particle i. The flowchart of 



Parameter i.e., velocity position initialization 



Velocity and position update with min and max limit 



Fitness value calculation with updated velocity and position 



Update J?„^ according to the fitness function 




Figure 1: Flowchart of SLPSO 

The particle i will learn (correcll its behavior vector) 
only if a randomly generated probability pi 
satisfies^ < J'i CO < < 1. In letail, ^^Lj(t + l>is 
generated as follows: 



fe) +ta(f).E.Ci.y(f) 



The swarm size m be determined as 



SLPSO algorithm is depicted in Fig 1. 



search dimensionality in the following form: 



(8) 



a function of the 



(9) 



where M is the base swarm size for the SL-PSO to work 
properly. 



IV. Results and Discussion 

The proposed problem of maximizing the DG 
penetration level is implemented in MATLAB and tested 
in IEEE 30 bus system. In this work, the population size is 
set to 20 and the maximum number of iterations is set to 
100. The single line diagram of IEEE 30 bus system is 
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depicted in Fig. 2. This network has six generators, 21 
loads and 41 branches. 

For IEEE 30 bus system the following assumptions are 
observed in this approach. 

> The boundaries of bus voltage limits are set at 0.95 
and 1.05 p.u., respectively. 

> The maximum numbers of DG units are taken as 
12. 

> The minimum and maximum values of DG units 
are limited between 0. 1 - 5 MW, respectively. 
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An optimal penetration level of 60.7950% is obtained 
for IEEE 30 bus system using SLPSO algorithm. The 
improved voltage profile after placing the DGs optimally 
is shown in Fig. 3. 
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Figure 3: Voltage profile improvement 



Consider the three phase fault is occurred at location 6, 



gh the bus is 9.6139 
used for protection 



Figure 2: Single line diagram of IEEE 30 bus system 

Base case load flow is performed and corresponding 
values of system power flows, bus voltages and line losses 
of test systems ar e observed. The bus corre sponding to tl ^ie 
minimum THD value is identified as optimal location of 
DGs. Hence, among all the calculated THD values the 
minimum 12 THD values are considered and their 
corresponding buses are identified as optimal locations for 
placing DGs. After finding these optimal locations, the 
efforts are focused towards the finding of optimal capacity 
of DGs, to be conn ected in these preferred locations, with 



then the fault current flowing throu 
p.u.. There are 6 relays are 
coordination, one is primary relay dnd other 5 relays are 
backup relays. If the fault is occurred at location 6, then 
the primary relay is used to detect the fault. If the primary 
relay is fails to operate then the backup relay is used to 
-detect-the fauit witlrirriiiimuiirq 
system from damage. The primary and backup relays 
opreting time is tabulated below in Table 2 and Fig. 4. 



of IEEE 30 bus system 



Table 2: Operating time of primary 



Operating times of relay in se : (p=primary, b=backup) 



Fault 



and backup relays 



the help of SLPSO algorithm. The optimal location and 
size of DG units are tabulated below in Table 1. 

Table 1: Optimal location and size of IEEE 30 bus system 



Location 



F6 



bl 



b2 



b3 



b4 



b5 



0.2600 0.7083 



1.3291 2.0120 2.7632 3.5896 



Total No 


Optimal 


DG sizes 


of 


Location 


in MW 


DGs 






1 


28 


4.9887 


2 


5 


4.7180 


3 


6 


4.9687 


4 


9 


4.9674 


5 


11 


4.7298 


6 


15 


4.4436 
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Figure 4: Primary and Backup relay operating timesi 

The comparison of real power losses without and with 
DG for IEEE 30 bus system is tabulated in Table 3. After 
placing the DGs, the real power loss is reduced from 
2.4438 MW to 1.6772 MW. In this case, the values of 
percentage reduction in real power losses are nearly 
31.3696%. 

Table 3: Comparison of losses for IEEE 30 bus system 



Losses 
without DG in 
MW 




Losses with 
DGs in MW 


Reduced 
Losses in % 


2.4438 




1.6772 


31.3696 



V. Conclusion 

This method of optimal placement and capacity of DGs 
using the SLPSO has been implemented for 30 bus systems 
The improper location and incorrect size of DG units can 
cause major disturbances in power system network. The 
DG penetration level is limited by fault current levels, 
protection coordination of relays, harmonic distortion, 
stability and power quality. 

The objective of the proposed optimization problem is 
to maximize the DG penetration level and loss 
minimization subject to power balance constraints and 
total harmonic distortion limits. The result shows that the 
maximum DG penetration level is limited by harmonic 
distortion level. Therefore, in this paper the bus 
corresponding to the minimum value of THD is identified 
as optimal location of DGs for the improved voltage 
profile and reduced losses in the systems. The proposed 
optimization method can serve as a good planning tool for 
the utility operator to optimally allocate DG of different 
types in order to achieve better penetration level. 
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